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Due to the advantages of a self-excited induction generator (SEIG), it plays a 
main role in sources of renewable energy, such as wind turbines (WT). The 
regulation of terminal voltage and frequency is poor under variable rotor 
speed and load conditions at stand-alone operation mode. The generator 
terminal voltage depends on the excitation capacitance which can be 
controlled by a capacitor bank and static voltage compensator. The 
dynamical model of the machine is described by differential equations in D- 
Q axes transformations of the synchronously rotating frame. Many models 
of analysis are proposed in the literature. In those models, several 
approximations are used to simplify the process of calculations, such as 
neglecting the iron core resistance, stray load resistance, stator and rotor 
leakage reactance, and magnetic saturation. In this work, a comprehensive 
dynamic model of the SEIG-WT is performed to analyze the system 
performance under transient and steady-state conditions. This dynamic 
model considers the effect of all machine parameters variation. New 
analytical formulas are used for to accurate calculation of minimum and 


maximum values of excitation capacitance and generator rotor cut-off and 
maximum speed. The dynamic model results are partially compared with 
experimental results, and accurate agreement is shown. 
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1. INTRODUCTION 

Due to the brushless construction of a self-excited induction generator (SEIG), it is also known as 
brushless asynchronous generator (BAG). This generator is a good candidate with the wind turbine for 
electric power generation, especially in an isolated and remote area as a stand-alone generation system with a 
self-excitation by a capacitors bank with the following advantages; brush-less construction, it can be excited 
without an extra power supply, easy to parallel operation, high reliability, and efficiency, no hunting, simple 
construction, and operation, it has a low unit cost, low running cost, low capital cost, ruggedness and 
robustness construction, better transient performance, absence of moving contacts, it has no DC power 
supply for field excitation, and low maintenance requirements. 

Although these are preferable features, induction generators have poor voltage and frequency 
adjusting on the generator terminal phases with the variation of load and speed. When the SEIG is directly 
connected to the national grid, it starts generating power when its shaft speed higher than the synchronous 
speed, which is limited by the grid frequency. The poor regulation of generator frequency may be adjusted by 
the speed-governor of the prime-mover. The blade pitch angle controller with the actuator adjusts shaft speed 
in the case of the SEIG-WT system. Generator terminal voltages can be regulated by excitation capacitances. 
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The stand-alone self-excited induction generator-wind turbine (SEIG-WT) system in an isolated 
mode of operation is shown in Figure 1. It consists of four components: wind turbine (prime-mover), SEIG, 
Excitation capacitor bank, and load. The reactive power supplied to the load is obtained by subtracting the 
reactive power consumed by the generator from the reactive power generated by the excitation capacitance, 
while the active load power equals the generated real power. 

Various dynamic models of SEIG analysis have been presented by researchers [1]-[7]. In all of these 
models, the cross-coupling saturation of magnetizing inductance, the dynamic saturation of leakage 
inductances, the influence of magnetic saturation on the iron core, friction, and windage resistance, and the 
stray load loss resistance are neglected, which leads to un-accurate results of those models. 

The present paper deals with dynamic modeling and simulation of stand-alone SEIG-WT that drove 
the system to feed an isolated area. The modeling of SEIG is synchronously rotating due to the voltages and 
currents being balanced. The dynamic magnetic saturation of magnetizing inductance, cross-coupling 
magnetic saturation of the magnetizing inductance, and the dynamic saturation of leakage inductances are 
considered in the dynamic model of the SEIG for the first time. Also, for the first time, the iron-core 
resistance and mechanical resistance are taken in the dynamic model in terms of air-gap magnetizing voltage 
by a polynomial curve fitting technique of the experimental setup measurements. The stray load loss 
resistance is considered in the model from measurements of the poly-phase induction machines [8]. Also, a 
new and simple method is used for the first time to calculate the minimum and maximum excitation 
capacitances, minimum and maximum speed and minimum and maximum slip of the SEIG based on the 
complete parameters of the generator, taking the effect of the saturated iron core and stray load loss 
resistances into consideration. 

The wind turbine dynamic model and simulation are adopted in this paper taking wind speed and 
blade pitch angle of the turbine as the main variables in the model. The results of the SEIG-WT simulation of 
the dynamic model are partially compared with experimental setup results to show the validity of the 
proposed model of the system. A variable speed, separately excited DC motor is emulated of the wind 
turbine, and the details of the parameters are given in the appendix. 

Finally, the paper is divided into six sections. Section one introduces the SEIG-WT system in the 
stand-alone and isolated mode of operation. Section two deals with the proposed equivalent circuit 
construction to derive the SEIG dynamical model including the standard tests to determine the equivalent 
circuit parameters of the SEIG for accurate modeling and analysis, taking the effect of saturation on the 
magnetizing inductance, leakage inductances, iron core resistance, and mechanical resistance. Also, this 
section deals with the effect of saturation on the cross-coupling magnetizing inductance and the effect of 
stray load loss resistance to perform a complete dynamic model of SEIG. Section three develops a simple and 
accurate analytical model for calculating the minimum and maximum excitation capacitance and generator 
rotor speed as well as slip at maximum torque. 

Section four deals with the wind-turbine dynamic model based on wind speed and turbine-blade 
pitch angle as the main variable parameters in the model. Finally, section five adopts the main results of the 
SEIG-WT model by comparison to these results with experimental results, and section six gives the main 


conclusions. 
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Figure 1. SEIG-WT system 
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2. DYNAMIC MODELING OF SEIG 

To complete the dynamic model of the SEIG, the model parameters must be measured or 
determined from the machine proposed equivalent circuit by performing the standard test on the machine, 
and the machine performance can be calculated at any load conditions without actually loading it. The effect 
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of machine inductance saturation on the iron core and mechanical losses are considered in the model. The 
effect of stray load loss on the machine efficiency is also considered in the dynamic model. The effect of 
cross-coupling magnetizing saturation on the machine performance is also considered in the dynamic model. 


2.1. Standard tests of the SEIG 

To model and analyze the SEIG, the parameters of its performance must be known at any load 
conditions, which are obtained from standard tests such as no-load, blocked rotor, DC resistance, and no-load 
at synchronous speed test. A setup approach has been adopted for the calculation of stator and rotor phase 
resistances, stray load loss resistance, unsaturated magnetizing and leakage inductances, unsaturated iron 
core resistance, unsaturated mechanical loss resistance, and the saturated magnetizing inductance, leakage 
inductances, iron core resistance, and mechanical (friction & windage) loss resistance. The dynamically 
saturated components of the equivalent circuit depend on the magnetic material type of stator and rotor iron 
cores. All the saturated inductances depend on the magnetizing current, while the saturated iron core 
resistance and mechanical resistance depend on the magnetizing air gap voltage. Figure 2 shows the complete 
equivalent-circuit components of SEIG, the equivalent-circuit parameters are given in the appendix. 
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Figure 2. The proposed SEIG equivalent circuit takes the effect of stray and mechanical losses 


2.2. Taking the effect of magnetizing and leakage inductance's dynamic saturation in the model 

The classical linear analytical model under the assumption of fixed or linear variation of the 
magnetic inductances leads to unrealistic results from the dynamic model due to the model failing to take the 
effect of saturation in both magnetizing and leakage inductances. Indeed, the nonlinearity due to the 
saturation effect is of crucial importance in the case of SEIG operation and performance, due to the narrow 
air gap of the machine, which leads to a non-negligible saturation effect. 

Therefore, for the high accuracy of simulation results, the iron core magnetic saturation and its 
nonlinearity have to be considered in the dynamic model simulation. In the model of SEIG, the dynamic 
saturation in magnetizing inductance and leakage inductances are the main factors for the stabilization and 
starting up of the generated terminal voltage. Hence, the inclusion of these inductance variations with 
magnetizing current and machine load current must be implemented in the dynamic model. 

Analytical expressions of magnetizing and leakage inductances from their practical relations with 
magnetizing and load currents, and a polynomial curve fitting technique are employed in Matlab, and the 
estimated functions are obtained as (1) to (3): 


Ly = 0.04765 * 15, — 0.5129 * I4 +2.18» I2, — 4.57 * I2, + 4.636 * Im + 0.25 (1) 


Le, = 4.824 * 1075 « 15, — 0.000714 * If, + 0.00414 « 13, — 0.01174 * I2, + 
0.01591,,, + 0.0156 Q) 


Ly, = 5.27 * 1075 « I5, — 0.0008987 * I5, + 0.00577 * I3, — 0.0176 * I2, + 
0.0251 * Ij, + 0.023 (3) 


by storing these functions in the Matlab Simulink of the dynamic model, the values of the saturated 
inductances can be updated at each step of integration. 


2.3. Magnetic saturation effect on the iron core and mechanical losses 
Nonlinearity in the iron core and mechanical resistances due to the saturation of the magnetic core 


can be taken into the SEIG model to increase the accuracy of simulation results. This condition is also 
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necessary for studying the voltage build-up in the isolated mode of SEIG operation since the linear model is 
not able of studying the performance of such systems. 

In general, the actual non-linearity in iron core and mechanical loss resistances can be found from 
the no-load test to assess the SEIG performance during the loading. Iron core resistance (Ric) and mechanical 
loss resistance (Ry) exhibit variations, but overall the iron-core resistance decrease with increasing air-gap 
magnetizing voltage while the mechanical loss resistance increases with increasing air gap voltage. Any 
change in load, speed, or excitation-capacitance will change the degree of saturation in the iron core, which 
consecutively will change values of Ric and Rs,, resulting in variable iron core and mechanical losses. These 
two resistances are determined experimentally at synchronous speed and modeled by the polynomial curve 
fitting as (4) and (5): 


Ric = —0.0001985 * V? + 0.111 * V? — 23.11 * V} + 840 (4) 
Rew = —6.167 * V? + 0.00419 * V? — 1.074 « V? + 127.4 * V, + 560 (5) 


the air gap magnetizing voltage (V,) can be calculated by adding the stator voltage drop in the stator 
impedance to the stator terminal voltage. 


2.4. Stray load loss resistance model 

To complete the dynamic model of SEIG for accurate performance calculations from the simulation 
results, all the losses determination must be considered in the model. This loss is represented by resistance 
(Rsee) inserted in the SEIG equivalent circuit as shown in Figure 2. This component is connected in series 
with the rotor resistance (R,-). Calculation of the stray load loss resistance requires performing a test at the 
rated full load, no-load, and the variable load of the machine. This loss is approximated by 1.896 of 
full load output machine power for small and medium-size machines and the stray load resistance can be 
calculated as [9]. 


Ree = 0.018 * Ry * (1 — Sm)/Sm (6) 


Where K, is the rotor resistance reflected on the stator side? S, is the full load slip at maximum 
torque and determined from the equivalent circuit components as will be shown later due to the minimum 
and maximum of SEIG excitation capacitance and speed depending on this slip. The physical representation 
and evaluation of stray load losses is a controversial and difficult problem in electrical machinery. IEEE 112- 
B Standard represents these losses by a parasitic resistance connected in series with rotor resistance in the 
rotor branch of the equivalent circuit, as shown in Figure 2. It follows that the stray loss depends on the 
square of the rotor current and the machine's full load slip. From the proposed equivalent circuit, the load 


resistance R,(1 — s) / s is derived as (7). 


R,—s)/s=1/ lu [3 -A4]- wal n 


7 f 


2.5. Cross-coupling magnetizing saturation 

The progress in the PWM technique, inverter-fed induction generator as voltage and frequency 
regulator as well as operation of an induction generator in isolated power system have led to an increase in 
the possibilities of including the mutual or magnetizing and leakage fluxes saturation effect in the dynamic 
model of SEIG. The saturation effect in the main flux of induction machines introduces a cross-coupling 
magnetizing saturation effect that is not predicted by the non-linear dependency magnetizing and leakage 
inductances modeling. The cross-coupling saturation is assumed to be entirely in the magnetizing inductance. 
Hence different models have been developed to obtain differential equations for analyzing the cross 
saturation effect [10]-[12]. In all of these models, the stator-constant reference frame model is utilized for 
analyzing the effect of cross-coupling magnetizing saturation, except in [10] the model cannot be solved due 
to it has six unknowns and four differential equations only. 

The curve of saturated magnetizing inductance is required to determine the static magnetizing 
inductance Lm and the dynamic magnetizing inductance (Ly). This curve is determined by the no-load test of 
the machine at synchronous speed for a variable supply voltage from 25% to 125%. At this test, the air gap 
voltage (Vo), the static magnetizing inductance (Lm), dynamic magnetizing inductance Lmy, and the cross- 
coupling magnetizing inductance (Lagq) are determined as [11]: 
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V; = Von — ls * (Rs + iXes) (8) 
Ly = fim) = = @ (9) 
Luy = FE = Ly + liml * E (10) 
Luaq = Luqa = ene, i (0) 


liml dim 


where Vpn is the machine terminal voltage (V). 

I, is the machine terminal current (A). 

R, = the stator resistance per phase (Q). 

Xp, = the stator leakage reactance per phase (Q). 

Am = the magnetizing flux linkage per phase. 

The magnetizing current (im) has real and imaginary components (ima and imq) and can be given as (12). 


im = ima +J img inl = ida + ta 


The components of the magnetizing inductance (Ly) are Lug and Luq and can be given as (13) and (14). 


Ls Lig dL 
= ima E Lq | alm 
bug = Lu + pag Mda = lM ttg (14) 
Also the components of the air gap voltage (V,) per phase can be derived as (15) to (17). 

dam dim dim à . 
Vja = oe = Luaq * ma + Luma * and = Ws * (Lua *lna + Luaq * Ing) (15) 

dàm dim dim : : 
Voq VUE Luaq * and + Lug * m = Ws * (Luq * ling + Luaq * ima) (16) 
y = [Và + Yea quU) 


Where Ama , A45 » Vga and Vj, are the D-Q axes components of magnetic flux linkage and magnetizing air 
gap voltage per phase respectively. w, is the angular frequency (rad/sec). 


2.6. Modified equivalent circuits of SEIG 

The iron core loss is normally modeled as equivalent resistance (R;,) that represents the power 
dissipated in the core identified practically from the no-load test. If the iron core loss is inserted in a dynamic 
model of the machine, a sixth-order model will be obtained. By appropriate modification of the parallel core 
resistance into serial equivalent resistance in the stator circuit, the dynamic model is reduced to the fourth- 
order model. The parallel branch can be modified into a series equivalent magnetizing inductance (Lm) and 
iron core resistance ( Ries). These series components of the parallel branch are calculated as (18) and (19). 


Ls Lu*Ric 
Lm = TRE (os Lo E) 
d 2 
Rics = .Ric*(0s LM)". (19) 


[Rit ws L?] 
Where w, = the voltage angular speed (rad/sec). 


The series iron core resistance in the parallel branch (Ries) is reflected in the stator side as a series 
voltage drop with an equivalent stator circuit, as shown in Figure 3, and can be calculated as (20). 
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This modified equivalent circuit of SEIG is a preferable model due to taking main and leakage magnetic flux 
saturation and iron core loss determination without increasing the model order. 
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Figure 3. Modified the proposed equivalent circuit of SEIG with reflection iron core resistance on the stator 
circuit 


2.7. Complete dynamic model of SEIG 

Modified equivalent circuit of SEIG in synchronously rotating reference frame as shown in Figure 
4. The 3-phase machine can be analyzed by an equivalent two-phase machine with Park's transformation [4], 
[10]. Both stator and rotor parameters are transformed into a synchronously rotating reference frame that 
moves with a rotating magnetic field. This model is very important and preferred for transient analysis and 
will be used in the present dynamic model of SEIG analysis by choosing the electrical machine currents as 
the main variables in the model. Figure 4(a) and Figure 4(b) show the D and Q axes equivalent circuit model 
of SEIG with excitation-capacitance, stray load loss resistance, and inductive load in a synchronously 
rotating frame of the stator and rotor parameters. 
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Figure 4. Modified equivalent circuit of SEIG in synchronously rotating reference frame, 
(a) D-axis model, and (b) Q-axis model 


From the model, the stator and rotor voltages at no load taking the cross-coupling magnetizing 
saturation are obtained as (21) to (30): 
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" : es ; dAma A 
Vsa = Rs *Usq + Les * + Rics * ima — Ws * Les * lsq + — Ws * Lmq * lng (21) 


Vsq Rs * lsq + Les * —À 4 Rics * ling + Ws * Les * İsa pe Um a4 Ws * Lima * ima (22) 


í lra 
EE — (ge * (Le, * — + Lm * imq) + 


— Wse * ud (23) 


Ura = 0 = (É, + Rsee) * Ga + Ley * 
dÀmd 
dt 


V =0= (R. + Rsee) * pq + Loy * T+ Dsp * (Le, * bra + Lina * ima) + 


En 
TUM + Wee * Ana (24) 
Zma = Lmaq * ua + Lima * us D- axis flux linkage (25) 
UM = Dmag * e Ling * (26) 
Lmaq = mi “ma * up cross-coupling inductance (27) 
ima = (isa + tra) D- axis current (28) 
imq = (sq + trq) Q- axis current (29) 
Wse = (Ws — wr) slip speed in rad/sec (30) 
Where w, is the generator shaft speed in electrical radian/sec. 
Lima and Lmg are D-and Q-axes components of magnetizing inductance. 
Let La = Lsa * Lya — Lina (31) 
jn weis. (32) 


where Lsa = Les + Lma D- axis stator self-inductance. 

Lra = Ley + Lima D-axis rotor self-inductance. 

Lsq = Les + Lmg Q-axis stator self-inductance. 

Lrq = Ler + Lmq Q-axis rotor self-inductance. 

Then substitute (25), (26), (31), and (32) into (21) to (24) and re-arrange these equations as a function of 
stator and rotor currents, the following four first-order differential equations can be obtained: 


disa _ Lra 


dt = La F [Vsa (Rs + Rics) * isa LI Rics T ing + Os(Lsq e E: + Lmq s irq) d 

L 

za xxm (Á, + Rsee)ira + Wse(Lrg x - + Ling ^ isq + Árqo)] — Tx Lmaq * Ws * img .. * (33) 
disq la, l ; 

dt = pe. * [vs = (Rs + Rics) i lsq PS Rics * lrq Ws(Lsa * Iq + Lma 3 tra) = 

q 

L A " E 5 L F 

ca [- (R, + Reee)trq E Wse(Lra * Íra + Lina * İsa + Árao)] = A * Lmaq * Ws * İma ...... (34) 
CUT NEA 


dt — ts * [7 (Rr + Rsee)íra + wselLra * frq + Ling * isq + Argo) — 
Vsa 


E * Ws * Í 35 
(Rs + Rics) * isa + Ws (Io i lsq t Lam * ing) — Rics * La La maq Ds Imq bd: ( ) 


eu 


- E [- (Á, + Rsee)írq — Wse(Lra * tra t Lina * İsa + Árao)] = 
Lmq 


E k Ling * Ws * ima ----(36) 


Vsq — | L 
Lq 


E") ke, + Rics) * lsq ni Ws(Lsq * İsa + Lma 2d tra) U Rics x irq » 
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The differential (33) to (36) are fed directly to the Matlab Simulink Block-set as a current space 
model to solve the first-order differential equations of the SEIG D-Q model without the need for matrix 
inversion in each step of integration as in conventional state-space models. The following auxiliary 
differential equations are used to complete the solution of the dynamic model of SEIG. The stator phase 
voltage components in the D and Q axes (Vsą and v,4) is calculated as a function of excitation-capacitance 
voltage in the D and Q axes (Veq and Veq) as (37) and (38): 


Vsa = Vea = Í (e + @, * Veq) dt + Veado (37) 


Vsq = Veq = Í (= — Ws * vea) dt + Veqo (38) 


where ig = iça — igg D-axis capacitance current component. 

icq = İsq — i94 Q-axis capacitance current component. 

igg and igg are the D-and-Q axes load current components calculated from the following differential 
equations in case of the resistive-inductive load: 


dipg 1 2 ° 

td — x = * * 

ae T 1, Uta 7 Re * iga) + ws * beg (39) 
Se = 1 4 (uy, — Ry * igg) -o*i (40) 
dt — m tq £ +q S td 


Where v;4 and Veq are the load voltage D-Q components and can be calculated as (41): 
Ved = Vea » Veq = Veq (41) 


R, is the resistive component of the load and L, is the inductive component of the load. 
Then the generated electro-magnetic torque inside the machine (T,) can be calculated as (42): 


T, = (2).(P,). Ga). (Ga a. Gata) (42) 


Where P, = the machine number of pole pairs. 
The angular frequency of the generated terminal voltage is to be calculated as (43): 


Qg — icq/(Ce * Vsa) (43) 
The D-Q components of the stator flux linkage are to be calculated as (44) and (45): 
Asa = SWsa —isq-R,).dt (44) 


Asq = S (Vsq — isq Rs). dt 
(45) 


The maximum stator phase voltage and current can be determined as (46) and (47): 


V, 


ph (m) = [v2 +v (46) 
Iph m) = fiz, + i (47) 


The active load power is given as (48): 
P, = 2 [vsa -iza + Voq Leq] (48) 
t~ 95 sd * td sq: "tq 
The total apparent power of the machine is calculated as (49): 
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S= 2 T Vpn(m) * Ion(m) (49) 


The load power factor can be calculated as (50): 


cos, = P; / S, (50) 


3. MINIMUM AND MAXIMUM SEIG SPEED AND EXCITATION VALUES IN STAND-ALONE 
OPERATION 

For calculation of the suitable value of the critical capacitor to initiate the SEIG-WT system 
successfully, many calculation methods are proposed in the literature [13]-[24]. 

The methods of excitation capacitors calculation, mentioned in the above literature are used 
inconvenient models with several approximations to simplify the steps of the calculation, such as a limited 
range of speed and resistive load, neglecting the iron core resistance, the stray load loss resistance, and the 
effect of magnetizing saturation. Also, most of the previous works use a numerical iteration method to 
determine the minimum and maximum excitation capacitance, but the iteration procedure takes a long time, 
and maybe a divergence occurs, and this process cannot be used online. 

The section contains a new, simple, and accurate analytical procedure that depends on the equivalent 
circuit parameters of the machine to determine the critical (minimum and maximum) excitation capacitor and 
critical (minimum and maximum) generator rotor speed values. The theoretical results from the proposed 
formulas are compared with experimental results, using the parameters of the induction generator configured 
in the laboratory to show the validity of the new method. The excitation-capacitance calculation is based on 
the resonance between the stator winding inductance and excitation capacitance at no load. Selecting the 
resonance frequency equals the machine base frequency at no load. The formula for calculation of excitation 
capacitance can be derived as (51): 


1 


Ce RETES (51) 


07 (Lit Les) 


In all literature, the authors have used this equation for determining the minimum excitation capacitance 
which is not a true assumption, the minimum and maximum excitation capacitance depend on the maximum 
and minimum generator rotor speed (w+) in electrical radian/sec. These minimum and maximum excitation 
capacitances are calculated as (52) and (53): 


t 


wre (Lm * Les) 


Ce (max) — (52) 
1 


wrm (Lm+ Les) 


Ce (min) = (53) 


Where w,, = the cut-off or minimum rotor speed in (electrical rad/sec). 
Wrm is the maximum (above synchronous speed) allowable machine speed in electrical radian/second 


Wrc = w; . (1 — Sm) (54) 
Wrm = Ws -(1 + Sm) (55) 


where w, = the angular-frequency in electrical. rad/sec. 
Sim = the slips at maximum torque. The slip is positive when the machine rotates under synchronous speed, 
and it is negative when the machine rotates above synchronous speed. 

The slip at maximum torque is calculated from the machine equivalent circuit analysis, taking the 
effect of iron core resistance (R;,) and stray load loss resistance (Rpp) as well as the dynamic saturation into 
account. By applying Thevenin's theorem on the proposed circuit in Figure 2, the following relations are 
obtained: 


R - [Ric XthmtRic ‘Rthm (RictRthm)| (56) 
mE [(Rict+Rthm)? +X tim] 


X ous Ric Xtnm (57) 
tht [tnm * Ric)? +X 
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Where R,;;=the total of the Thevenin's equivalent resistance of the circuit. 
Xtne=the total of Thevenin's equivalent reactance. 
Renm=the Thevenin's equivalent circuit resistance and is determined as (58): 


Rs X74 
[R2 TX es Xm)? ] 


Rinm = (58) 


Xtnm is the Thevenin's equivalent reactance taking the effect of magnetizing inductance only, and is 
calculated as (59): 


Xen, — Xml Ri Xs Xes+Xm)] 
thm [R2 Gto * X0)? | 


(59) 


The maximum torque of the machine occurs at maximum power transfer between the stator and rotor, and 
this condition will occur when the rotor resistance component (R,/s) equal to the other components of the 
circuit as (60): 


Ky : 
T = Rent Rege + i(Xene + Xor) (60) 


S 
The slip at maximum torque is obtained as (61): 
Sm = Rr / (Rene + Rsee)® + Xene + Xer)? ]^ (61) 


For approximation results, if the magnetizing inductance and iron core resistance are neglected the slip at 
maximum torque becomes (62): 
Sm = R, / [(Rs + Ry)? + (Xes + Xy)? ]? (62) 


4. THE TURBINE-GENERATOR DYNAMIC MODEL 

The wind turbine is modeled by relations of power coefficient (C;) as a function of both tip speed 
ratio (A) and blade pitch angle (£). The value of / is fixed at its maximum values, and the power coefficient 
(C5) will become maximum also. The tip speed ratio is the ratio of angular turbine rotor speed to linear wind 
speed at the tip of turbine blades [25], [26]. It can be given as (63): 


à= o, .R/V, (63) 


R is the blade length. V,, is the speed of wind in (m/ s) and w+ is the shaft speed in mechanical rad/second. 
The swept wind power in Watt is calculated as (64): 


B, — 2.ps. A Va (64) 
Where pa is the air-specific density in (Kg / m?). 
A is the swept area by the blades in (m?). 
The turbine output power (P;) can be calculated as (65): 

P, — C, * P, (65) 
The turbine torque in (N.m) can be written as (66): 

T, =C,* B, (66) 
Where C, is the torque coefficient and given as (67): 


WS EU (67) 


A suitable way to control the power obtained by the wind turbine is to adjust the blade pitch angle 
(B). The blade pitch is similar to the gate valve in hydroelectric power turbines, except that its response is 
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much faster. At a downwind speed (V, € 13 m / s), the value of f is fixed at approximately Zero degrees 
and not required to be adjusted, while the turbine operates at maximum power. If the wind speed becomes 
greater than 13 m / s ), the turbine power will highly increase due to it being increased with cubic wind 
speed, and it must be adjusted by controlling the blade pitch angle to prevent the turbine from destruction. 
The controller will feed the required or demand p (fq) to the actuator, and the actuator (stepper 
motor) will adjust the required pitch angle to reduce the turbine speed and power at fixed rate values. For a 
realistic actuator response to the pitch angle, the actuator uses the time constant integrator, so that the pitch 
angle is from 0° to 20° with a rate of change +1° per second. The actuator can be modeled as (68) [27]: 


B=— S Ga- Pm) dt (68) 


Where Pm is the measured pitch angle by the actuator and 7g is the actuator time constant in deg/sec? 
The power coefficient (Cp) in terms of tip speed ratio (A) and the turbine blade pitch angle (£) in degree is 
obtained as (69): 


(69) 


Cy (B) = C, * 5 — C3 -£ = c, * NS + Cs A, Ai = | 0.0848 t +A+B*+0.08+8 +A | 


B3-2.8*10-3*B—-0.035*A+1 


where C4 = 0.5176 , C; = 116, C3 = 0.4, C4, =5, Cs = 21 and C6 = 0.0068 
The generator shaft speed in (electrical. rad/sec) can be calculated by integrating the first-order differential 
equation of the motion as (70): 


= Px (Ta — Te — D * Wm) (70) 


Where wpm is the generator shaft speed in (mechanical rad/sec) and can be calculated as (71): 
Om = 0. [ P, (71) 


where J is the total inertia of the turbine-generator system in (Kg. m?), and calculated as (72): 


] 7 (&*1,) (72) 


where J; is wind turbine moment of inertia and Jg is the generator moment of inertia, and can be obtained 
from the deceleration (retardation) test. G is the gear-box ratio and can be obtained as (73): 


G=% = 


wt Tsh 


(73) 


Tsn is the torque on the generator rotor ( load torque ) in N.m developed by the turbine and can be calculated 
as (74): 


Tt 


T =Ë (74) 
Where T; is the aerodynamic torque of the turbine shaft and can be calculated as (75): 
P 
T, — t/a, (75) 


D is the viscous friction torque coefficient in [N.m/ (rad./sec.)] and can be calculated as (76) [28]-[30]: 
D = Pg / we, (76) 


Where: Pr, is the mechanical power loss. 
T, is the generator electromagnetic torque calculated by equation (42)? 


5. RESULTS AND DISCUSSION 
In the experimental setup, a 3-phase, squirrel cage, 220V, A- connected induction machine is joined 
to a wind turbine and acts as a prime-mover. Test results and simulation results in Matlab/Simulink of the 
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SEIG-WT system are presented at no load initiation of the self-excitation. The SEIG is driven at 1525 r.p.m 
(above synchronous speed). Then a connected delta excitation capacitance bank of 40 uF was switched 
across the SEIG stator terminals. The induction machine and wind turbine specifications and parameters are 
given as: 

SEIG specifications and equivalent circuit parameters: 

3-phase, squirrel cage rotor, 220 V, 50 H; , A- Connected with the following parameters: 


Pout = 1500 W ,1, = 6.6 A, N, = 1500 r.p.m , cos = 0.8, P, = 2, J, = 0.03 K,.m? 
D = 0.003 N.m/(Rad./sec.) , Cex(min) = 30 Hp, Rs = 5.10, R, = 3.5 Q, Cut-off speed = 1200 r. p.m 


Design class B, Insulation class-F 
Wind turbine specifications: 


Pout(max) = 1750 W, Blades = 3, pg = 1.225 g, Cp(max) = 0.4, R= Im (radius), G = 1.25 (gear ratio). 


Variations in magnetizing inductance and stator and rotor leakage inductances as the main 
influencing parameters have been considered into account by the experimental non-linear relationships 
between magnetizing-inductance and current and between leakage inductances and load current as in (1) to 
(3). Variations in iron core and mechanical resistances are taken into account by experimental non-linear 
relations between these resistances and air gap magnetizing voltage as given in (4) and (5). Building- 
generator terminal voltage as shown in Figure 5. Variation of permissible cut-off and maximum rotor speed 
has been tested experimentally, as a parameter to the generated voltage and calculated analytically by (54) 
and (55). Figure 5(a) and Figure 5(b) show the experimental and MATLAB simulation results for starting- 
generator terminal voltage under transient conditions. 


Figure 5. Building-generator terminal voltage, (a) experimental result and (b) MATLAB simulation result 


Building-generator terminal voltage as shown in Figure 6. From these figures, it can be shown that 
very good agreement between the experimental results and theoretical simulation due to the accurate 
simulation model, which considers into account all the SEIG-WT system parameters and variables. 
Figure 6(a) and Figure 6(b) show the experimental and Matlab simulation results of the steady-state terminal 
voltage of the generator. Figure 7 shows the MATLAB result of the electromagnetic torque of the generator 
at starting-up conditions. 
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Figure 7. Matlab simulation results of electromagnetic torque of SEIG during the starting-up condition 


6. CONCLUSION 

In this dynamic model of the SEIG-WT system, all important parameters of SEIG are considered 
without introducing extra differential equations to the model, due to the modification in the equivalent circuit 
of SEIG. The dynamic saturation of magnetizing-inductance, cross-coupling magnetizing-inductance, and 
leakage-inductances are taken in the model due to these parameters are the main variables that influence the 
SEIG performance in steady-state as well as transient conditions. The influence of magnetizing saturation of 
the air-gap voltage on the iron core loss and mechanical loss, as well as the influence of stray load loss, are 
considered in the dynamic model for accurate prediction of the system performance. New, simple, and 
accurate analytical formulas are derived for the calculation of minimum and maximum excitation 
capacitances and rotor speed of SEIG. These formulas are obtained from the comprehensive equivalent 
circuit of the machine and based on maximum power transfer theory. 
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